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Abstract 

A new Langevin equation with a field-dependent kernel is proposed to deal with bot- 
tomless systems within the framework of the stochastic quantization of Parisi and Wu. The 
corresponding Fokker-Planck equation is shown to be a diffusion-type equation and is solved 
analytically. An interesting connection between the solution with the ordinary Feynman 
measure, which in this case is not normalizable, is clarified. 



1 



One of the interesting and appealing applications of the stochastic quantization (SQ) of 
Parisi and Wu [1] is found in its application to systems characterized by actions unbounded 
from below (bottomless systems). An attempt to stabilize and eventually to quantize bot- 
tomless systems is of great interest in several important fields of physics, such as Euclidean 
Einstein gravity. To deal with bottomless systems on the basis of SQ, Greensite and Halpern 
[2] made use of the Fokker-Planck equation w.r.t. a fictitious time t 

^P[0;t] = i/[0]P[0;t], (la) 

with an action S unbounded from below. In this case, the distribution functional 
(Feynman measure) cannot be an equilibrium solution nor even a stationary solution, for it 
is not normalizable. Because the negative semi-definiteness of the Fokker-Planck Hamiltonian 
H can be proved irrespectively of the boundedness of the action, they proposed the interesting 
idea that a possible way to stabilize and to quantize such bottomless systems in Euclidean 
space-time is such that the probability distribTition in the Euclidean path integral formulation 
is given by the lowest (normalizable) cigcnstatc of H. This corresponds to using the lowest 
eigenstate of H of the form 6^*^°^', instead of e^'^, as an effective distribution functional 
in evaluating expectation values. It is also shown that such an effective distribution can 
indeed be extracted by fixing both initial and final configurations [3]. This procedure avoids 
runaway solutions which arc due to the unboundedness of the action S of the Langevin 
equation. Notice that their proposal for such bottomless systems corresponds to a departure 
from the traditional approach that the weight functional in the path integral formulation be 
given by e~^, with S a classical action. Remember also that their Langevin equation does 
not have an equilibrium distribution because the naive candidate e"*^ is not normalizable 
and hence does not belong to the spectrum of the Fokker-Planck Hamiltonian H: Every 
eigenstate belongs to its negative-definite eigenvalue. 

Recently Tanaka et al. proposed another interesting way of stochastically quantizing 
bottomless systems [4]. Unlike the above attempts [2,3], they have pursued the possibility of 
reproducing the desired probability distribution e"*^ in the equilibrium limit of the fictitious 
stochastic process even for bottomless systems. Their starting point is the following Langevin 
equation 

with a ficld-dcpcndent positive kernel A'[0]. Here ?7 is a Gaussian white noise with the 
statistical properties 



{r]{x, t)) = 0, {r]{x, t)r]{x', t')) = 2S% - x')5{t - t') 

2 



(3) 



and we adopt, here and in what follows, the Ito-interpretation [5] for stochastic equations. 
It is straightforward to derive the corresponding Fokker-Planck equation from the above 
Langevin equation (2), and the Fokker-Planck Hamiltonian is then given by 

It is clear from this equation that for any positive definite kernel functional the spec- 
trum of H is still negative semi-definite so that the system relaxes to a unique equilibrium 
state whose distribution is given by e^*^, if it is normalizable and if the spectrum has a 
discrete zero. If the classical action S is not bounded from below, the distribution e"*^ is 
not normalizable, and can never be either an equilibrium solution or a stationary solution of 
the above equation. Nevertheless, a better choice of the kernel K may result in the situation 
where the stochastic variable is effectively confined to a finite range, so that the distribu- 
tion can be regarded as a solution of (4) normalizable within such a range. Their choice 
of the kernel K, for example, is 

K[(f)] = e-^4[0] (5) 
when the action S is given by the difference of two positive monomials 

S2[(p] = \j d'^x<p{x)[-d'^ + m2]0(a;), 54[0] = d^x<p\x), 
This choice results in the following form of the Langevin equation 

which has a damping drift term. They have argued [4] that this drift term effectively confines 
to a finite range so that the desired distribution e"*^ can be derived in the equilibrium limit 
of (7). Their numerical simulation of the above kerneled Langevin equation for a simple zero- 
dimensional model (wrong-sign x'^ model) could indeed reproduce the distribution e"'^'-*'* = 

1 2 2 A 4 

g-^m X +jx ^-^ equilibrium [4]. 

Seeking a possible way of reproducing e"*^ in the equilibrium limit for bottomless sys- 
tems is challenging and very appealing. However, their proof is crucially dependent on the 
plausible argument that the effective range of the stochastic variable may be finite, which 
makes the apparently nonnormalizable distribution normalizable, but the soundness of 
which is only verified numerically. 

In this short paper, I propose another way of effectively reproducing e~^, starting from 
the above Langevin equation (2) with a slightly different kernel K. It will be shown that the 
resultant Fokker-Planck equation can be solved exactly at any fictitious time t < oo. It turns 



(6a) 

(A > 0). (6b) 
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out to be a diffusion-type equation, i.e. the spectrum of the Hamiltonian H is continuous and 
no equihbrium hmit exists. This might be considered an apparent drawback, but it wih also 
be shown that the asymptotic form of the probabihty distribution is well approximated by 
the desired form normalized in the finite range of the stochastic variable, corresponding 
to a finite fictitious time. These results, though of course not contradictory to it, form an 
interesting contrast with the previous analysis [4], where a discrete spectrum of H is implicitly 
assumed and the equilibrium distribution is numerically found to be well approximated by 

Throughout this paper, we exclusively consider zero-dimensional potential models whose 
classical action S{x) is unbounded from below, for notational simplicity. 

Let us start our discussion with the simplest possible model, i.e. the wrong-sign Gaussian 
model 

S{x) = -\mx^, (m > 0), (8) 

and show that an exact solution to the Fokker-Planck equation (la) with the Hamiltonian 
(4) is obtainable if we choose an appropriate kernel K{x). First we divide the above action 
into two positive monomials as in (6a), 

S{x)^Sq{x)-S,{x), (9) 

with 

Sq{x) — ^mx^, Si{x) — mx"^. (10) 
Then following Ref . [4] , we are naturally led to a kernel 

X(x) = e-^'(^) =e-"^^'. (11) 



Though the above separation ((9) and (10)) seems quite arbitrary, it will be shown that this 
choice is crucial for the ensuing Fokker-Planck equation to be solvable and is just an example 
of a general scheme K — e^^ which in this case is equal to (11) because 25" = —S^^ (see (29) 
and (32) below). The exphcit forms of the Langevin and the Fokker-Planck equations are 

i ^ -e-^i(^)5^(x) + e-2^'(^)77 = -e-^^'mx + e'^^^^, (12) 

and 

P(x; t) = ^e-<^) (A + P(.; t) = |e-»^-' (| - m.) P(.; t), (13) 

respectively. As we know that the corresponding Fokker-Planck Hamiltonian has negative 
semi-definite eigenvalues, we put 



P(x;i) = e-'^V^(^)p(x), 
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(14) 



with a real parameter c. The above Fokker- Planck equation (13) is then reduced to an 
ordinary differential equation for P{x) 

-c2e"^^'p(x) =(^-mx] ^Pix). (15) 
\dx J ax 

Now we change the variable from x to f{x) in the above equation and obtain 

|(/(a:))'^ + {f"{x) - mxf\x))^ + c^e-^-'l P{x) = 0. (16) 

For the equation to be solvable, we require 

f'{x) oc e"'^'/2 and f"{x) - mxf{x) = 0. (17) 
It is clear that this condition is identically satisfied for f'(x) oc e^^ or 

f{x)^a Tdye^y'/^ + P (18) 



with two real parameters a and (3. Notice that the transformation between x and the above 
/ is one to one and the domain of / is again [— oo,cxd]. This choice of f{x) drastically 
simplifies the equation (16) and makes it solvable. In fact, P now satisfies 

|«'^ + c2|p(x)=0, (19) 

which is easily solved to give 

P{x) = Re{e''Io'ye"''"^') or Im(e-r^^^'"''''). (20) 

At this stage it is important to notice that the normalization condition, which is the only 
condition that the probability distribution P is subject to, 

dxP{x]t)= dxe-''^e~^^'-'>P{x) = 1 (21) 

-00 ^ — 00 

seems to be satisfied for any real c 7^ because at large \x\, which corresponds to large 
1/1, rapid oscillations of P (see (20)) may make the integral convergent despite the ever- 
increasing factor e~^^^\ No condition has been imposed on the value of c, which implies 
that the system has a continuous spectrum. This will be verified shortly. 

Thus we have arrived at the conclusion that the solution of the Fokker-Planck equation 
(13) is in general of the form 



J— 00 '- ^—00 



(22) 



where the function A{c) in the second equahty has been re-defined to absorb the constant 
term in f{x). Specification of A{c) certainly corresponds to a specific choice of the initial 
distribution P{x; 0). It is not difficult to see that the simplest choice of A{c) does correspond 
to the most fundamental initial condition. In fact, if we choose a constant A — l/27r, we 
have, at the initial time t = 0, 

P{x;0) = r dce"'^(*)Re 

In the last equality we have made use of the one-to-one correspondence between x and 
^Qctye^y and the monotonically increasing character of the latter w.r.t. x. 
In this way, we finally obtain the solution of the Fokker-Planck equation (13) 

/oo 2 , n/ \ 

dee-" V^(^)Re 
-oo 

subject to the initial condition P(x; 0) = 5{x). It is easily confirmed that this solution is 
normalizable, i.e. that J^^dx P{x;t) = 1. As a trivial generalization, the solution subject 
to the initial condition P{x; 0) = Pq{x) appears to be 

P{x-t)= da;o^=e-'^(^-)e ^V-o''^' ^Po{xo). (25) 
i-oo yAnt 

It is instructive to observe that this analytical solution obtained for the simplest case can be 
expressed in terms of S only 

P{x; t) = r dxQ e-^(^)e~^te '^^'^"^^''O Pq{xo). (26) 

J-oo V47ri 

This suggests that there may in general exist an analytical solution. This is indeed the case, 
as will be shown explicitly later. 

It is remarkable that we can obtain analytical solutions to the Fokker-Planck equation 
(13), corresponding to the nonhnear Langevin equation (12). It should be stressed again that 
the solution decays at large t and therefore no equilibrium limit exists and its behaviour is 
thus quite different from that expected in the ordinary case: It is rather similar to that of 
diffusion processes, which again implies that the spectrum is continuous. If, for example, we 
choose a trivial kernel i^' = 1 in (13), no nonlinearity appears and the standard technique is 
applicable even though the potential has the wrong sign [6]. In this case, the solution Pk=i 
subject to Pk=i{x; 0) = 6{x — xq) is 

P^=i(a;;t) = y^^^^^— ^32^exp 



27r 



dy e 



d{x). (23) 



27r 



-S{X) 



(24) 



m{xe — xq) 



2(1 



-2mt 



-mt 



(27) 
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whose asymptotic behaviour at large t shows the discreteness of the spectrum 



Pk=i{x; t) ^-mxl/2^-mt_ (28) 

V ZTT 

Now let us turn our attention to more general cases. No specific form of the classical 
action S is assumed here, except for its unboundedness from below. The problem is to find 
an appropriate kernel K which makes the Fokker-Planck equation solvable, as in the previous 
example. If we write a positive kernel K as 

K{x) = e-^(^) (29) 

using a real function W{x), the function P already introduced in (14) satisfies 

|(/(^))'^ + (fix) - [W'{x) + S'{x)]f{x))^ + c2e^(-)| P(a:) = 0. (30) 

In this case, solvability requires that 

/'(x) oc e^(^V2 and f{x)-[W'{x) + S'{x)]f'{x)^0. (31) 

These conditions are clearly satisfied by the following / and W 

W{x) = -2S{x), (32) 
fix) = a' r dy e-^^y^ + (3'. (33) 

J 

Having found the appropriate kernel and the transformation function, we are therefore able 
to solve the Fokker-Planck equation, which in this case is written 

P(a;;t) = ^e2^w(^ + 5'(a;))p(a;;t), (34) 

just as in the previous example. The solution P{x;t), satisfying the initial condition 
P{x; 0) = Po{x), is straightforwardly found to be 

/•oo 1 or ^ r dve-^^yA'^ 

P{x-t)= dxo^=e-^^'=h ^'V^o^y^ ^ Po{xo), (35) 
i-oo yAnt 

as expected from (26). 

Finally, let us discuss possible implications of these results in numerical simulations. One 
usually keeps track of the stochastic variable, whose time development is governed by the 
Langevin equation 

i; = e2^(^)5'(x)+ 6^(^)77 (36) 



when the kernel is given hy K — e^^ . Let the total action 5" be given by the difference 
between two terms 

S{x)^Sq{x)-S,{x), (37) 

where is a positive monomial of the highest order, assumed even, and Sq{x) the remaining 
terms. The above Langevin equation is approximated for large x by 

i ^ -e-25'(^)5;(a;) + e'^'^^^Ty, (38) 

and for small x by 

X ~ e2^o(^')5^(a;) + e^o^^^r;. (39) 

Thus the kernel K = e^^ supplies a restoring force for large x, but a diverging force for small 
X. To understand the behaviour of the stochastic variable x{t) starting from a finite value 
a;(0) = xq, and the approximate form of the probability distribution P, consider only the 
finite range \x{t) — xq\ < -R. This range may be viewed as an actual domain of the variable 
in numerical simulations; (no infinite range can be realized in practical simulations.) Next 
we define tn. as a time scale satisfying 

Then at t = the second exponential factor of the exact solution (35) is expressed as 

,->fe(/;*«-''")' = e-te^'^'OVte-n"*--*')' (41) 

and is considered to be of the order of unity for any x in the range \x — .tqI < R. Under 
this estimation, the probability distribution at time is well approximated by a naive 
distribution normalized in this range 

P(a;;ta)~^— , Zr= dye-'^y\ (42) 

This implies the possibility in the actual simulation of the Langevin equation (36) of repro- 
ducing the desired probability distribution at a certain finite time tj^, normalized in a 
finite domain determined by the details of the simulation. Further investigation is necessary 
to draw more definite conclusions. 

The author acknowledges useful and informative discussions with Profs. I. Ohba, Y. Ya- 
manaka and K. Okano and Drs. S. Tanaka, M. Mizutani and M. Kanenaga. 
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